Long range orders in multiferroics 
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We proposed that in multiferroics there exists a third long range order besides the electric polar- 
ization and magnetic order. This long range order reduces the symmetry of the spatial part of the 
wave functions of electrons. Thus the cancellation in the "spin current" model can be avoided. As 
a result, the expectation value of electric polarization will be larger by an order of magnitude. We 
have derived a new form of electric polarization P ~ — Q x (sj x Sj+i) where Q is the wave vector 
of this long range order. 

PACS numbers: 75.85.+t, 77.80.-e, 75.25.Dk 



Multiferroic materials [ij have several coexisting long 
range orders. Kimura et. al. found very interest- 
ing properties in TbMnOa. The ferroelectric transition 
temperature coincides with a magnetic transition which 
was identified later to be the transition from sinusoidal 
magnetic order to spiral order Q. It was also found Q 
that with sufficiently strong magnetic field, the electric 
polarization switches direction. These are clear experi- 
mental evidences that the electric polarization and mag- 
netic orders are closely related. Electromagnetic field is 
one of the fundamental interactions and electricity and 
magnetism are unified. It should not have been surpris- 
ing that electric and magnetic properties are coupled. 
However, in this case there seems to be no charge cur- 
rent nor time-dependent electric field, nor varying mag- 
netic fiux in the system. So the fundamental question is 
by what mechanism electricity and magnetism are cou- 
pled. As the list of multiferroic material grows rapidly, 
their mechanisms and the varieties of physical proper- 
ties also multiply Van der Brink and Khomskii 
classified them into two types Q- Type-I multiferroics 
usually have high ferroelectric temperature but weak cou- 
pling between electric polarization and magnetic order. 
In type-II multiferroics, there are strong electric polariza- 
tion and magnetic order coupling. With a few exceptions 
such as Ca3Co2-a;Mna;06 [s] and the scenario proposed 
by Sergienko et. al., Q, these compounds are originated 
from spin-orbit interaction and possess spiral spin con- 
figuration. This is the focus of this work. Katsura et. al. 
(KNB) proposed the spin current model of multifer- 
roics in which the coupling between electric polarization 
and magnetic order is due to spin-orbit interaction. The 
mechanism can be interpreted as the coupling between 
electric field and spin current, or the Aharonov-Casher 
(AC) effect ll| and it gives polarization 



transfer energy, e'l and 62 are the unit vectors of two 
adjacent spins and ei2 is the direction of bond, el = 
e 1/ 'ip*^{r)r^p{r)d?'r\, is the expectation value of electric 
polarization of the hybridized state of d-orbitals of transi- 
tion element atoms and p-orbitals of oxygen atoms. The 
factor ei x 62 fits nicely with spiral spin configuration. It 
predicts that the electric polarization is perpendicular to 
the screw direction. Furthermore, it was able to explain 
the results of Kimura et. al. applied magnetic field 
can flip the direction of spiral spins and hence, that of 
electric polarization. The generalization to a bulk system 
with active Cg orbitals was also performed • 

Though the spin current model has been successful in 
many aspects, it has the quantitative problem. For ex- 
ample, assuming V « O.SeT^ and A k, 2eV, I has to be of 
the order O.lnm in order that P sa lOOOfiC/m^, the mag- 
nitude in TbMnOa. This raised doubt of the relevance of 
the model [13|, . First-principle calculation by Mala- 



P ~ -el 



ei2 X (ei X 62) 



(1) 



where V is the hybridization energy, A is the charge 



shevich and Vanderbilt [15[ indicated that in TbMnOa 
the electric polarization of purely electronic nature is 
an order of magnitude smaller than experimental data. 
This problem has to be addressed before the spin current 
model can be viewed as an accomplished model. This is 
also the subject of this work. 

The electric polarization in Eq.([T|) which has the factor 
{V/A) to the third power, is the remains of catastrophic 
cancellation which comes from the symmetry of the sys- 
tem. In KNB's words: "the dominant term comes from 
the difference of the normalization factors between two 
perturbed states." The electric polarization exists only in 
the systems without inversion symmetry. In many com- 
pounds it is the spiral spin configuration which breaks the 
symmetry. However, as one can see from the calculation 
of KNB, this effect, to the leading order, is only manifest 
in the spin part of the wave functions and it bears the 
factor ei x 62. The spatial part of the wave functions, 
to the leading order, still possess space inversion symme- 
try. Hence, the expectation values of electric polarization 
suffers cancellation. There are other evidences that the 
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lower the symmetry of the system, the greater the electric 
polarization. In Q , it was found that atom displacement 
enhanced electric polarization in E-type antiferromag- 
netism. The results in ref. provide even stronger sup- 
port. Their calculation showed that the spin-orbit inter- 
action is indispensable to the existence of electric polar- 
ization while lattice distortion enhances electric polariza- 
tion by at least ten-fold. However, Bridges et. al. [3] did 
not find evidence of lattice distortion within experimen- 
tal accuracy, which is 5 x 10~'^A. It must be said that a 
displacement of ions approximately 1 x 10~^Ais suf- 
ficient to give the electric polarization found in TbMnOs. 
It will be even smaller for the displacement of Cu^"^ ions 
(4 X 10~^A) in the compounds with smaller polarization, 
like LiCu202 [17]. At this stage, experimental data can 
neither support nor rule out lattice distortion contribu- 
tion in TbMnOa and LiCu2 02. Therefore, it is likely that 
in the compounds with large electric polarization such as 
TbMnOa and CuO there is an extra factor which 
reduce the symmetry. As a result, the spatial part of the 
system also has lower symmetry to avoid the catastrophic 
cancellation. 

In this work we proposed that there should be a third 
long range order besides the electric polarization and 
magnetic order. It can be structure distortion, charge or- 
der, or orbital order (00). 00 occurs in transition metal 
compounds with crystal field, Jahn- Teller or GdFeOa- 
type distortions [l9| . For example in manganites, ip^^^-r^ 
and '03y2_r2 orbitals are favored on two inter-penetrating 
sublattices due to crystal field [l^l- However, other ef- 
fects such as superexchange interaction and spin-orbit 
interaction can also change the orbital configuration. In 
the spin-orbital model [l9| different combinations of or- 
bitals can give different exchange energies. Hence, there 
are various competing mechanisms. As a result, the 00 
should not have such a simple structure. We can pro- 
vide another argument for the existence of 00. In 3d- 
transition elements the spin-orbit coupling strength is of 
the order O.OSeV or 600K. It is much higher than the Neel 
temperature. This means that the spin and orbital angu- 
lar momenta are aligned in opposite direction long before 
the magnetic order is established. If there is experimen- 
tal evidence of say, sinusoidal magnetic order, it is also 
an indication of 00 with a wave vector commensurate 
with that of spins. 

Even if in a system without 00, structure distortion, 
such as cooperative Jahn- Teller effect or GdFeOa-type 
distortions can also play the role of the third long order. 
The d-orbitals are anisotropic. The magnitude of oxygen 
bond angle can influence hybridization. The third possi- 
ble candidate is charge order or its multipole expansion. 
For apparent reason it can cause an additional periodic 
potential and hence, affects wave functions. 

Here is a brief summary of this work. By introducing 
a third long range order into the system, we show that 
the electric polarization can be enhanced by an order of 



magnitude. A new form of electric polarization is derived. 
This long range order can reduce the spatial symmetry of 
the spatial wave functions. Therefore, the catastrophic 
cancellation is avoided. More specifically, long range or- 
ders with wave vector Q affect the hybridization energy. 
The state of crystal momentum k mixes with the state 
of k±Q. In evaluation of electric polarization one has to 
sum over entire bands up to the Fermi surface for multi- 
ferroic insulators. The cancellation will not be complete 
if the wave function is a mixture. 

To see the origin of the catastrophic cancellation we 
start with an U(l)xSU(2) gauge-invariant action which 
will give the equation of motion as the Pauli equation: 

/t2 
d^xiihroo^p - — Y.iDki^TiDk^') 
k=l 

~ jtc t ^.^^n (2) 

with x'^ — ct and the covariant derivatives Dq = d / dx^ — 
{ie/hc)Ao + {ie/2mc^)J2aPoaCT'' and Dk = d/dx'' + 
{ie/hc)Ak + {ie / Amc^) J2 PkaO'"' are introduced in order 
to show clearly the charge and spin currents. A term of 
EM field is also present. Here, we use the notations of 
Frohlich and Studer [21] with a little simplification. The 
SU(2) fields poa = and pka = J2b ^kabE^ are related 
the components of magnetic and electric field. However, 
it should be noted that the action is correct only to the 
order 0{m~^). Above action illustrates the origin of AC 
effect in which the electric field is coupled to moving spins 
or magnetic moments. By taking variation of pka^ mul- 
tiplying by Ekab and summing over a, and assuming the 
system is charge neutral, we obtain 

Eb -he rr,h d e ^ , „ , 

^ = ^^ZJ , A[{-jrT + -^kWekaba''4' 
Aire om'^c'^^ — 'a,k=i % cix'" c 

+ (efca6^"V')*(-/^ + -Afe)V} (3) 

I ox'" c 

where a term of the order m^^ is discarded. The r.h.s. is 
clearly the spin current. It is the gist of the spin current 
model of multiferroics. From the equation D ~ E + AttP 
one finds that when there is no free charge, the electric 
polarization P = E/Air can be generated by spin current. 
Nevertheless, we can see the reason why the electric po- 
larization given by the spin current model usually is small 
in magnitude. In insulators such as multiferroics, there 
is no movement of free charge carriers and thus no spin 
current on the r.h.s. of Eq.ljSj. KNB had to invoke su- 
perexchange interaction. This is exactly the reason why 
there is catastrophic cancellation because only by going 
to the higher orders can one find spin current. The can- 
cellation can be avoided if there is a long range order in 
the system which can give the wave function modulation, 
produce significant derivative of ^ and hence, give rise to 
large the electric polarization. 
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Our model is quite simple. A crystal of transition ele- 
ment atoms with valence d-orbitals and oxygen atoms 
with p-orbitals are considered. The transition metal 
atoms have spiral spin configuration and 00 so that the 
j-th transition metal atom has the combined spin and 
orbital state 



cos{q-Rj/2) 

[A + Bcos{Q-Rj)]i^di{r-Ri. 
Bs\n{Q-Rj)4,d2{r-Rj) 



(4) 



where ipdi{2){'f' ~ Rj) is the first and second d-orbital of 
the transition element atom at site q and Q stand 
for the wave vectors of spin order and orbital order re- 
spectively and is the angle between x-axis and the spin. 
We assume that there are two active orbitals like those 
of Bg orbitals and have reserved the freedom of 00 with 
non-orthogonal orbitals, like \?>x^ — r^) and |3j/^ — r^). 
B is approximately a constant. If on the other hand, the 
third long-range order is lattice distortion and there is no 
00, then A = 1 and B ~ 0, (or one can set Q — 0.). The 
effect of lattice distortion will appear in the Hamiltonian 
in next paragraph in the form of modulated hybridization 
energy. 

The hybridization energy is affected by the orbital 
states and lattice distortion due to the anisotropics of 
the orbitals. It is also influenced by spins. The spins 
of the hopping electrons are aligned with localized spins. 
Hence, the Hamiltonian has the form 

+ ^[AVi+ BVcosiQ • Rj - ai)\ cos(Q' • Rj) 

n.n. 

[e""''/^ cos0Jc4■Cp^,^t + e"^^^ sin^jC^Cp*,;;] + H.c. 
+ A^/;-s,- (5) 

where I is the index of the oxygen atoms (position Rj-\-fi) 
in the basis, and 6j — q ■ Rj/2. Vi and V/ are the hy- 
bridization energies of the p-orbital of the l-th oxygen 
atoms with the first and second d-orbital of the transi- 
tion metal atoms, V = \A;+Vp^ and cos a; = Vi/V. 
The last term is the spin-orbit interaction. The factor 
AVi + BV cos{Q ■ Rj —ai) in hybridization comes entirely 
from 00. The modulation BV cos{Q ■ Rj — a/) contains 
contribution from the first and second orbitals with the 
factors cos a/ and sin a; respectively. If there is no 00, 
then A ^ 1, B = 0, ai = 0. The factor cos(Q' • Rj) 
comes from the displacements of atoms. For example, 
if the oxygen atom between two transition metal atoms 
move laterally such that the bond angle of TM-O-TM 
is reduced, then the hybridization will be reduced for 
eg electrons. Hence the hybridization energy acquires 
the modulation of lattice distortion (wave vector Q'). If 



both 00 and lattice distortion are present, then our case 
is viable only if they are commensurate to each other. 
In that case, two effects are superimposed on each other 
and the hybridization energy will has the form above. In 
the following calculation we will consider the case with 
00 only. The extension to additional lattice distortion 
effect is not difficult. 

The Hamiltonian without spin-orbit interaction can be 
diagonalized by making the substitution: Ri i — Ri + fi, 
PUA = e-**/2 cos(g-i?,^,/2)cp,,/t + e*^/2 sin(g • i?,,i/2)cp,,^ 



and pii^i = e"^/2cos(g • Ri,i/2)cpiji - e 



ict>/2 



sm{q 



Rij/2)cpi^i^, and taking Fourier transform. For our pur- 
pose, it suffices to assume that Sd — £p > \V\ and analyze 
the perturbed wave function. The resulting wave func- 
tion can be expressed with the following tight-binding 
wave functions: "0^ ^(r) — ''^^ e'''''^^'^j and i^pij. {'r) = 

'^e'^-^'il)p{r- Rj ~ n)(e~''^/2 cos 6'^ , e'-^/^ sin6'j)'^and it 



is 



E. 



B^ 



p,Lk+CQ 



Op L 

(6) 



where ipp{f — Rj — n) is the p-orbitals of the oxygen 
atoms at site Rj + fi, = ±1, and C is a normalization 
constant. 

We apply the spin-orbit interaction to the wave func- 
tions in Eq.Q as a perturbation and gain a term 
— {\/A)l ■ s'i'j where A = Jh — Ac/ with Jh being the 
Hund's coupling energy and Ac/ the crystal field splitting 
between e^, and t2g orbitals. (The perturbed states have 
minority spin and t2g orbitals.) Now we are in position 
to calculate the electric polarization with Eq.([6]): 



P 



c.c. 



(7) 



The unperturbed wave functions do not contribute to 
electric polarization because of cancellation. It is due 
to crystal structure inversion symmetry despite of spiral 
spin configuration. The first term vanishes after sum- 
ming over k. It is again due to the inversion symmetry 
of the spatial part of the wave function. The net electric 
polarization comes from the second term in Eq.([7]). The 
reason the catastrophic cancellation is avoided is that 
with the third long range order, the energy bands open 
gaps at 2k- {Q+G) = (Q+G)^ where G is a reciprocal lat- 
tice vector. The gap is comparable to the hybridization 
energy if 00 or lattice distortion is significant. Differ- 
ent mini-bands not only have different energies but also 
have different linear combinations of wave functions. In 
other words, the symmetry of the entire band is lower. 
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FIG. 1: (color online) Schematic graph showing the direc- 
tions of spins (red arrows), electric polarization (grey arrow) 
q (black arrow) and Q (blue arrow), q and Q can be parallel 
to each other. See text. 



Introducing phonons into the system [l5| has the same 
effect. The importance of Q is manifest. Calculation of 
the contribution of the highest mini-band gives 



P ! 



^el'ABW cos a 



7rA(£d - Ep) 
[cos(Qa;ao) sin(Qa;ao/2) sin(ga;ao/2) cos(/) 
-I- cos((5a;ao) sin((5yao/2) sin((7j,ao/2) sin^]. (8) 

where ao is the lattice constant and I' — J ipl^ (f — 
Rj)z^x{r — Rj — fi)d^r. We define cosa = | cosazj, note 
that for — y^ orbital T4 — —Vy and for 3z^ — r^ orbital, 
Vx = Vy, so that cosa^: 
limit where ao ^ 0, we found 



ibcosaj,. In the continuum 



2erVX 



Q X (sj X Sj+i). 



(9) 



Here, s" is the unit vector of the spin. In a crystal with a 
long-range order other than spins and electric polariza- 
tion, the wave vector of the long-range order Q, emerges 
in the expression of electric polarization in place of the 
bond direction in Eq.([T]), as shown in Fig. 1. 

Since we are concerned with the magnitude of the elec- 
tric polarization, we now make a practical estimation: 
\Vi\ ^ O.beV, Ed - £p w A leV, A « O.OSeF and 
AB « 0.5. Note that above calculation is purely of hy- 
bridization. It does not involve double occupancy in the 
transition metal atoms and hence, the on-site Coulomb 
repulsion does not have any effect. For this reason Sd — Sp 
and A are smaller than those calculated with LDA-I-U 
density functional theory. The magnitude of electric po- 
larization is |P| ~ 3000 - 5000(/7 A)/^C/m2. Since the 
order of magnitude of /' is a fraction of A. The result of 
our calculation is compatible with the experimental data 
of say, TbMnOg. 

On reflection of Eq.(I3|), one sees clearly the origin of 
the electric polarization in Eq.([3]). The ip in Eq.(I3]) is just 
that in Eq.([n]). The scattering of electrons with the third 
long range order, which can be seen from our calculations 
in Eq.((31ini), mixes different plane waves and hence, en- 
able ip to sustain current. With the magnetic moments 



provided by spiral spins in leading order, the spin current 
on the r.h.s. of Eq.(|31) gives rise to electric field on the 
l.h.s. which in turn, causes electric polarization in the 
form — With linear response formalism 



|i?-^|n)(n|7^|0) 

Eo ~ En 



(10) 



one finds that in Eq.(|ni) and Ea. ([TU|) has the same 
form. We have the following correspondence: En — Eq ^ 
A, ^ • 7^ gives spin-orbit interaction with E in Eq.([2]) 
and the spin operator in Eq. gives the now well-quoted 
factor X Sj+i- Since spins and orbitals are closely re- 
lated, wave vectors of magnetic orders and orbital order 
very likely have the relation Q = rq where r is a rational 
number. Our calculation showed that only incommensu- 
rate 00 can suppress catastrophic cancellation. If not, 
00 can be accounted for with a enlarged unit cell, and 
the electric polarization still suffers catastrophic cancella- 
tion. This can explain why magnetoelectric phenomenon 
usually occurs with incommensurate magnetic orders. 

Our result has several implications. For example, it 
was found that [3 SI] there is modulated electric polar- 
ization (antiferroelectricity) in multiferroics. Its absolute 
magnitude is much greater than that of the net polariza- 
tion. The existence of 00 makes things more complex. 
If the modulated polarization and 00 are commensurate 
to each other, then an enhancement can occur. A related 
phenomenon is electromagnons [23| in multiferroic mate- 
rial 24 1 . It was suggested [1^ that the dominant term 
is of the form P — J^^ij^i ' ^j- According to Moriya 
(23| . ~ I' J/ A where the superexchange interaction 
J {V in Moriya's notation) is of the order {V'^/A)^/U 
with U being the on-site Coulomb repulsion. Comparing 
to this P, our polarization in Eq.ijHl) more favorable 
or the spin-orbit coupling (0.05eV) is relatively stronger 
compared to exchange interaction. The additional long 
range order can give rise to extra magnon excitations. 
These are of purely electronic origin. The wave vector of 
the magnons can also be different. It is equal to Q for 
one-magnon processes and k and k±Q for bi-magnon pro- 
cesses. Its effect will be manifest in optical spectra, sus- 
ceptibility and electric polarizability measurement. This 
will be related to the vector spin chirality [lij Sj x Sj+i. 
The study on this subject is under way. 

In conclusion we have proposed that with third long 
range order in the system, the catastrophic cancellation 
of electric polarization can be avoided. The resulting 
magnitude is comparable to experimental data. The long 
range order can be either 00 or cooperative Jahn- Teller 
distortion, but it has to be incommensurate. Lastly but 
perhaps most importantly, the cause of electric polariza- 
tion and multiferroics is firmly established and it is spin 
current coming from the spin-orbit interaction. 
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